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Abstract. We discuss the the GLV Reaction Operator formalism for com- 
puting the non-abelian energy loss of jets propagating through hot and dense 
Q\ ' matter. We incorporate our results on induced gluon bremsstrahlung together 

with the effects of nuclear shadowing and multiple scattering in a two com- 
' ponent soft+hard description of heavy ion reactions at RHIC energies. We 

demonstrate that good agreement between data and theory can be achieved in 
| the measured moderate pr < 5 GeV window at ^/s = 130 AGeV and we extend 

. our predictions for the high px part of the hadronic spectra. We focus on the 

perspectives of using jet tomographic methods for probing not only the density 
and geometry of the quark-gluon plasma but also its parton composition at 
various center of mass energies. 
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1. Introduction 

Coherent scattering effects on induced photon/gluon emission have been studied in 
both abelian and non-abelian gauge theories. Destructive interference of photon 
radiation from initially on-shell fast electrons propagating and multiply interacting 
in matter was first discussed by Landau and Pomeranchuk [ [jj and independently 
by Migdal [ |2). Photon bremsstrahlung off asymptotically prepared jet states was 
found to be suppressed in the infrared (oj — > 0) region relative to the naive Bethe- 
Heitler limit [ [| of radiation resulting from independent collisions. Those results 
have been experimentally confirmed at SLAC [ || . 

The non-abelian generalization of Bethe-Heitler radiation was discussed by 
Bertsch and Gunion [ |s| in the presence of a single scattering center. In a QED- 
like scenario the case of multiple interactions was first investigated by Gyulassy 
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and Wang ||. The non-abelian nature of QCD coherence was later employed by 
BDMPS [ 0j to show that in the high energy regime the mean energy loss of asymp- 
totic jets is proportional to the square of the length of the medium, i.e. AE ~ L 2 . 
Those results were obtained for "thick" plasmas with a very large number of soft 
scatterings. Discussion of gluon bremsstrahlung in the path integral formalism can 
be found in Refs. [|, |. 

The GLV reaction operator approach [ [nj to non-abelian parton energy loss 
relies on a systematic expansion of induced gluon radiation associated with jet pro- 
duction in a dense QCD plasma in terms of correlations between multiple scattering 
centers. Analytic expressions for the induced inclusive gluon transverse momentum 
and light-cone momentum distributions are derived to all orders in powers of the 
opacity of the medium, \ = No g jA = L/X g . The analytic solution to all orders in 
opacity generalizes previous continuum results by allowing for arbitrary correlated 
nuclear geometry and evolving screening scales as well as the inclusion of finite kine- 
matic constraints. In comparison to data we use numerical results corrected up to 
third order in \ that allow us to extend jet quenching computations down to parton 
energies E ~ 5 GeV. For gaining theoretical insight, however, it is useful to resort 
to analytic formulas that neglect some of the kinematic constraints but capture the 
essential features of parton energy loss. The dominant first order energy loss can 
be written as 



where Cr is the second casimir in the representation of the jet and the /x 2 (r)/A(r) 
is the local "transport coefficient". In a hot and dense medium in local thermal 
equilibrium /i 2 (r) = ATta s T 2 (r) and perturbatively one can express /^ 2 (r)/A(r) = 
{CtCaIAj^^-kcc 1 p(r), where CV is the color charge of the target and is the 
dimension of the adjoint representation. This formulation allows to account for 
the realistic dynamical time evolution of the systems created in ultra-relativistic 
nucleus- nucleus collisions. 

Jet tomography is the QCD analog of conventional X-ray tomography that 
exploits the attenuation of high energy jets produced in nuclear matter [ [llfl . We 
use the term "jet tomography" , inspired by a seminar by Istvan Lovas [ |l2j] entitled 
"Vector meson tomography of the QGP" , to contrast the particular advantages 
of high pt PQCD probes of the density evolution of partonic matter formed in 
A+A reactions. The dependence of the non-abelian energy loss on the density 
and geometry of the medium as seen in Eq. ([j]) provides a rigorous theoretical 
framework for the jet tomographic analysis discussed in Refs. [ O, [l^, [H], pr| . 
The experimental discovery of a factor of ~ 3 suppression of moderate pr ~ 4 GeV 
7r°'s in central Au + Au reactions by PHENIX [ |l8| and the discovery of large 
transverse asymmetries in non-central collisions for p^ < 5 GeV by STAR [ have 
confirmed that the high pt frontier at RHIC does in fact provide a wide range of 
new physics opportunities. 

We here present the results of our calculation of the suppression of high pt 




(1) 



Jet Energy Loss in Hot and Dense Matter 



3 



particle spectra relative to the binary collision scaled pp reference baseline and the 
Pt and centrality dependence of the p/ir ratios [ We find good agreement 

between data and theory for jet energy loss driven by initial gluon rapidity density 
dN 9 jdy ~ 800. One of the descriptions of soft parton production, based on the 
saturation models [ |2^] , suggests slow variation of the number densities of quanta per 
unit rapidity achieved in the early stages of A+A reactions with s/s and therefore 
a slow variation of the quenching factor. In contrast, at SPS energies moderate pt 
neutral pions were found to be enhanced by a factor ~ 2 [ |2l[ . We propose that 
one possible solution to this puzzle may be related to differences in the composition 
of the non-abelian plasma (quarks vs. gluons), an important detail that has so far 
been ignored in jet tomographic studies. 



2. Induced Gluon Radiation 

The GLV formalism [ [h]] for multiple elastic and inelastic interactions inside nu- 
clear matter expands the differential probability of observing a final state jet or 
jet+gluon system (described by a set of quantum numbers {a}) in orders of the 
correlations between multiple scattering centers, i.e. P({a}) = ^2 n PW({a}). In 
the high energy eikonal approximation a major simplification occurs as a result of 
the well defined path ordering of sequential interactions inside the medium. It is 
therefore possible to build all relevant classes of amplitudes by subsequent insertion 
of single Born or "direct" (D) and double Born or "virtual" (V) interactions. At the 
probability level the recursion relation is generated by the reaction operator [ |Ic| ] 
R = D*D + V* + V. The virtual corrections V< + V to the naive elastic scattering 
component &D ensure unitarity in the GLV formalism. The approach described 
here is quite general, i.e. within the framework of the approximations stated for an 
arbitrary initial condition described by an amplitude A , pW({a}) cx Aq (R) n A . 
The technical part of the calculation includes solving for the color and kinematic 
structure of the the "direct" and "virtual" operators. Recently this was illustrated 
via a computation of the elastic broadening of jets propagating in nuclear matter [ 

For the case of jets produced at finite time to in heavy ion reactions the solution 
for the double differential gluon radiation intensity was given in Ref. [ [l0| 



~2 C(l,.»,n) ' B( m +x,—,n)(m,—,n) x 
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where C^...^) = |Vkk>g(k-qi - ••• -q„) 2 , TS( m +i,-,n)(m,-,n) = C( m+1 ...„) - 
C( TOl ..., n ) are the Bertsch-Gunion terms [ [|, and W(fc,..., n ) = 1/(2x15 IC^...^) | 2 ) are 
the interference phases that produce the non-abelian analog of the LPM effect [0,0]. 
The destructive interference pattern can be seen at a differential level in Eq. (0). 
In the calculation of the gluon bremsstrahlung we use a normalized color-screened 
Yukawa potential leading to v 2 (q^) — /i 2 /(7r(q 2 + /J 2 ) 2 ). One unexpected result was 
that the series in Eq. (j^) converges very fast [ [H], E3]. Even more surprisingly, 
the first term was found to give the dominant contribution and by integrating over 
the transverse momentum k and light-cone momentum fraction x of the gluon one 
recovers the result of Eq. (Q) . In computing the spectrum of gluons and the energy 
loss of jets we take into account numerically corrections up to third order in the 
opacity \ and the evolution of the density of the system due to the longitudinal 
Bjorken expansion. 



3. Hadronic Spectra at Moderate to High pr 

The computed energy loss spectrum is employed in the calculation of the quenched 
spectrum of hadrons. A jet of flavor c and transverse momentum p c produced in 
a hard PQCD scattering a + b ^ c + d is attenuated prior to hadronization by 
the radiative energy loss to p* = p c (l — e). This shifts the hadronic fragmentation 
fraction z c = p h /p c to z* = z c /(l - e). 

The invariant distribution of hadrons reduced by the energy loss in A + A 
collision at impact parameter b is then given by 

(IN AA f _ 

Eh ~(Pp~ = Taa ^ X! J dx i dx 2 d 2 k T , a d 2 k T ^g(k T , a )g(k T ^) 

abed 

fa/A{Xi,Q )fb/A{x 2 ,Q ) — t — / deP(e,p c ) ' , (3) 



dt 

where Taa{^>) is the Glauber profile density. The fragmentation functions D h / C {z ) Q 2 ) 
is taken from BKK [ |24|] and the structure functions for f a /A( x iQ 2 ) ( we take the 
GRV94 LO [ ^H) include the isospin dependence. Nuclear shadowing, fcx broaden- 
ing and Cronin effect are taken into account as in [ ^6], . The modification of 
the fragmentation functions due to gluon radiation can be seen in Eq. (|J) and is also 
discussed in [ ^9). The FFs carry the information of the medium induced jet energy 
loss but the observed suppression of high px is a result of the full calculation that 
folds in a variety of nuclear effects as well as the shape of the initial jet spectrum. 

In order to extend our calculations to lower transverse momenta of hadrons 
we discuss a phenomenological soft component motivated by the string and baryon 
junction picture. We patrametrize here this component as follows: 



dyj^~ ~dy~ { } 2^(T«(b)) 2 



„„ W • (4) 

7T,K,p, 
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Such seemingly "thermal" distribution of hadrons is observed already in e + e~ col- 
lisions and may be related to Gaussian fluctuations in the string tension [ p0[ . 
As in string models the soft component is assumed to scale with the number of 
participants (N part ). In Eq. (Q) we also account for the possibly different mean 
inverse slopes T a for baryons and mesons. In the junction picture [ |3l| , the 
large T B may arise from the predicted [ smaller junction trajectory slope 
a 1 j ~ a 'n/3- This implies that the effective string tension is three times higher 
than l/(2ira' R ) « 1 GeV/fm leading in the massless limit to (p\) j — ?>(p\)b.- In 
terms of the string model the factor three enhancement of the mean square pt is due 
to the random walk in p^ arising from the decay of the three strings attached to the 
junction. Naively, we would thus expect T B ~ \/2>T' K predicting an approximately 

constant J (p^) B for all baryon species. 
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Fig. 1. The ratio of charged hadron and ir° multiplicities to the binary collision 
scaled pp result is shown from [ [L8| . The curves utilize the GLV quenched hard 
spectrum and the string and baryon junction soft component Eq. (Q) . 

Fig. 1 shows the ratio Raa of the differential particle yields relative to the 
binary collision scaled pp for inclusive charged hadrons and neutral pions. We ob- 
tain reasonable agreement with data if the jet energy loss is driven by initial gluon 
rapidity density dN 9 /dy ~ 800. This number is lower by ~ 30% from existing 
estimates. The difference in the suppression factor of tt° and 0.5(h + + h~) is un- 
derstood through the possibly different baryon and meson production mechanisms 
in the moderate high 2 < pt < 5 GeV window. In our calculation pion production 
becomes PQCD dominated for p? > 2 GeV and correspondingly suppressed by the 
jet energy loss. In contrast, baryon production in the region of interest is dominated 
through the junction mechanism by baryon transport in rapidity and moderate pr- 
This accounts for the different suppression of neutral pions and inclusive charged 
hadrons as seen in Fig. 1. At large px baryon production also becomes perturbative, 
leading to a common suppression factor. This is manifest in the pr > 5 GeV region 
of Fig. 1. 
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Fig. 2. The centrality dependence of p/n is predicted for three different central- 
izes. Solid (dashed) lines correspond to A 1 (A 4 ^ 3 ) scaling of the baryon junction 
component. The ratio of p and ir~ fits to PHENIX data on central reactions is 
shown for comparison. A boosted thermal source (dashed line) is also shown. The 
left/right panel reflects negative/positive hadrons. 

The interplay between soft and hard physics at RHIC is possibly most clearly 
seen in the px differential baryon/meson (Rb) ratios. Our predictions for the cen- 
trality and the pr dependence of p/tt~ and p/ir + are given in Fig. 2. An extended 
discussion can be found in Refs. [ [H], [53|. In central collisions the interplay be- 
tween the anomalous baryon component and the quenched PQCD component of 
pions leads to maximum of Rb near px ~ 3 — 4 GeV/c. At large pt > 5 — 6 GeV/c 
we predict a gradual decrease of Rb below unity consistent with the the PQCD 
baseline calculations [ |l7|. In Fig. 1 we have included through error bands the 
factor of ~ 3 uncertainty in the fragmentation functions into p, p at high p^. The 
solid and dashed curves reflect the difference between the N part and N^^t scaling 
of the junction component. 

In peripheral reactions the size of the interaction region as well as the initial 
density of the medium decrease, leading to a reduction of energy loss. The absence of 
quenching reduces the observability of the anomalous component and the p/ir ratio 
may stay below unity for all pr- The case of peripheral reactions is hence similar to 
pp collisions. The experimentally testable prediction of the model is therefore that 
the maximum of the Rb = p/ft ratio decreases with increasing impact parameter, 
decreasing participant number, or equivalently decreasing dN ch /dy. The reduction 
of Rb at large pr is also an important prediction of our computation and is not 
seen by other models attempting fits only in regions of px where experimental data 
already exists. 
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4. Jet Tomography of the QGP Partonic Composition 

Two extreme scenarios for the composition of the quark-gluon plasma have been 
extensively discussed in the literature. In one scenario, based on saturation models [ 
, the hot and dense matter created in A+A collisions is interpreted as purely 
gluonic degrees of freedom. Alternatively, quark coalescence models [ |3f| advocate 
the creation of dense quark-antiquark matter as a result of the vary fast decay 
of massive gluons into quark degrees of freedom [ ^6). We here propose that jet 
tomography is a powerful tool that can distinguish between different scenarios. 
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0.1 



10 100 1000 

s 1 ' 2 [AGeV] 




Fig. 3. The effective fractional energy loss for E = 6 GeV jet is plotted versus 
yfs for two extreme scenarios. 

Fig. 3 illustrates the idea behind the QGP composition tomography. We con- 
trast the case of gluon jets propagating through gluon plasma to the case of quarks 
propagating through quark-antiquark plasma at center of mass energies ranging 
from y/s = 20 AGeV at SPS to y/s = 200 AGeV (the top RHIC energy). The 
variation in the effective fractional energy loss (using an approximately constant 
density renormalization factor Z ~ 0.4 — 0.5 [|23|) is small ~ 50% due to the small 
change of the initial density of quanta per unit rapidity (this can be most easily 
seen from the naive picture of parton-hadron duality). In contrast, the changing 
plasma composition and the jet representation lead to a factor of (Ca/Cf) 2 — 5 
difference. While neither of the extreme scenarios may hold for the broad range 
of collision energies in question, it is possible that the partonic composition of the 
system created in heavy ion reactions evolves as a function of <Js, leading to pro- 
nounced differences in AE e ff(y/s)/E as indicated in Fig. 3. Faster than expected 
increase of energy loss with T/T c was also discussed in a Polyakov Loop model [ 

The quenching of neutral pions is shown in Fig. 4. In this particular example 
we study only the effect of energy loss and we have not included nuclear shadowing 
and Cronin effect. We have also not included a phenomenological soft component 
to mock the part of phase space that is not perturbatively accessible and is possibly 
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sensitive to higher twist effects. In Fig. 4 the medium that drives the quenching of 
high px jets at RHIC yfs = 130 AGeV is reinterpreted as composed of 60% gluons 
and 40% quarks and antiquarks. At yfs = 200 AGeV Au+Au the number of quanta 
per unit rapidity increases only by ~ 15% [ [57j. The dashed line illustrates the 
suppression factor Raa for this slightly increased rapidity density of soft partons 
when unchanged QGP composition is assumed. We note that the effect is even 
smaller than naively expected (~ 15%) due to finite kinematics effects in computing 
AE/E and hadronization. Detectable increase ~ 30% in the 7r° suppression factor 
is only possible if one assumes a transition to almost completely gluon dominated 
soft background as shown in Fig. 4. 

At smaller center of mass energies (y/s = 20 AGeV) we have studied the case 
of jets propagating through quark- antiquark plasma. We have used phenomenolog- 
ically a smaller mean energy loss renormalization coefficient Z = 0.2 to account for 
the possibly bigger energy loss fluctuations (more detailed studies of Z are needed) . 
Those effects tend to strongly reduce the effective energy loss. The observed suppres- 
sion factor Raa, however, convolutes the apparent decrease of AE e ff/E with the 
much steeper pt differential hadron distributions (^). As a result a still significant 
suppression of neutral pions (Raa ~ 0.5) remains, which is inconsistent with the 
current SPS results [ |l| . It is important to improve the jet tomographic calculation 
presented in Fig. 4 to take into account Cronin effect and nuclear (anti) shadowing 
in a more quantitative comparison to the WA98 data as well as to investigate the 
cause for the unusually small energy loss at SPS. 
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Fig. 4. The 7r° suppression factor resulting from the non-abelian energy loss of 
jets is shown. 



5. Conclusions 



We have shown that in a two component soft+hard model that incorporates jet 
energy loss computed in the GLV formalism [ [l(| one can account for the suppression 
of the inclusive charged hadrons and neutral pions seen in the y/s = 130 AGeV RHIC 
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run. Jet tomography is shown to be an important probe of the partonic composition 
of the hot and dense medium created in the initial stages of heavy ion reactions. 
In particular, at CERN SPS energy of y/s ~ 20 AGeV large energy loss induced by 
gluon medium on gluon jets seems excluded by data. 
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Notes 



Recent reanalysis of the WA98 data (M.M. Aggarwal et ai, Eur. Phys. J. C23, 
225 (2002).) suggests that the enhancement estimates should be reduced by 
~ 30%. 

All recent approaches [ |?], |[ ^|, that discuss the LPM effect in QED and 
QCD employ the eikonal approximation together with a model of well separated 
scattering centers /iA ^> 1. For the case of A+A reaction this condition sets a 
natural upper limit on the opacity \ of a medium of finite small size L ~ 5 f m 
and is suggestive of the dominant role of low order correlations between multiple 



scatterings as demonstrated by GLV [ 10 
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